Vortex lattices: from the hydrodynamic to the quantum hall regime 
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We describe the excitations of a vortex lattice in rapidly rotating trapped Bose-Einstein Con- 
densates, focusing on their evolution from the hydrodynamic to the quantum Hall regimes. We 
observe the dramatic changes and show that single particle states are an increasingly dominant sig- 
nature of the spectrum. We interpret the behavior by analogy with the integral quantum Hall effect. 
Furthermore, we compare our results to recent experiments at J1LA and find excellent agreement. 
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The angular momentum of a rapidly rotating trapped 
Bose-Einstein condensate (BEC) is carried in singly 
quantized vortex lines. At equilibrium, the vortices typ- 
ically form a triangular array which rotates with the 
condensate Q, Q] - the Abrikosov lattice. Recent ex- 
perimental and theoretical efforts to understand the dy- 
namics of this vortex lattice have focused on the calcu- 
lation 0, H S B B 13 and measurement H @ of two 
particular branches of the spectrum commonly described 
in terms of the associated Tkachenko and inertial modes. 
The Tkachenko waves which are governed by the small 
but finite shear modulus (C2), mainly involve fluctua- 
tions in the velocity held of the condensate and appear 
to an observer as elliptically polarized transverse waves 
moving through the lattice (Fig. QJi) . On the other hand 
inertial modes which mainly involve density fluctuations 
similar to hydrodynamic shape oscillations are charac- 
terized by the bulk modulus (Ci). The finite size of the 
vortices ensures that both modes are intimately coupled 
and they can be shown to occur as two spectral branches 
of the coupled fluid-lattice equations in a hydrodynamic 
treatment 0. 

In the BEC, the angular frequency of the trap rotation 
(f2z) plays the role analogous to a uniform magnetic field 
applied to a type-II superconductor or to a field effect 
transistor (FET) as in the integer quantum Hall effect 
(IQHE) When investigating vortices in the trapped 
BEC, it is natural to ask whether the connection to the 
IQHE can provide insight and lead to new results. In 
this letter we concern ourselves with the transformation 
of the excitation spectrum as we transition from the stiff 
Thomas-Fermi (TF) regime to the mean field quantum 
Hall regime (QHR). We highlight edge modes associated 
with the boundary of the condensate which appear as 
a ring of vortices precessing the apparently undisturbed 




FIG. I: Azimuthally symmetric (3,0) a) Tkachenko mode 
and b) inertial mode, c) Surface wave 



ground state (FigHJ;) . These modes are particle-like exci- 
tations that have an intimate connection to Landau levels 
in the IQHE and are shown to be especially important 
for the BEC as the QHR is approached. Our identi- 
fication evokes an interpretation of the Tkachenko and 
inertial modes as analogues of extended and bound sin- 
gle particle states fundamental to the IQHE. Finally, we 
use this framework to discuss lattice meting and vortex 
nucleation. 

One exciting feature of trapped BECs is the wide range 
of rotational speeds over which vortex lattices may be 
produced. As f2 increases from to the radial trapping 
freqency u± , centrifugal effects cause the particle density 
to change and we encounter three physical regimes [3( 
which may be identified by three length scales defin- 
ing the lattice. These are the radius of the condensate 
in the x — y plane (R), the healing length(f) and the 
intervortex spacing (e) defined by £ = 1/ \/8ima and 
e 2 = 27rh/(\^3mQ). Here a and n represent the s-wave 
scattering length of the trapped species and density of 
particles at the trap center, e and £ are the coherence 
lengths of Tkachenko and inertial modes and togther with 
the wave vector k they determine whether an excitation 
has a collective or particle-like character. In the stiff 
(e 2 » £R) TF regime (£ <§C e <§C R), which occurs at slow 
rotation, the time scale of density fluctuations is much 
shorter than that of rotation and the system responds 
to perturbations as a rigid body. At faster speeds we 
encounter the soft (e 2 <C £i?) TF regime indicating that 
while the mean field interaction continues to dominate 
the properties of the lattice, the dynamics of the den- 
sity had become more sluggish and thus important in 
the response of the lattice to perturbations. At ultra- 
fast rotation we enter the QHR (e < £ < R) , here the 
chemical potential fi is dominated by the cyclotron en- 
ergy. In the QHR, the lattice is described by the primary 
occupation of the almost degenerate lowest Landau level 
(LLL) and is expected to undergo quantum phase 
transitions into a myriad of strongly correlated vortex 
liquid states 01 • The change in character of the excita- 
tions as a function of the relative sizes of the length scales 
(£, e, R) is an important theme in the transformation of 
the spectrum and the subject of this paper. 
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FIG. 2: a) The dispersion curves of azimuthally symmetric Tkachenko and inertial waves for lattices in different physical regimes. 
For each plot the scales labelled u>i and lot on left and right ordinate are for inertial and Tkachenko modes, respectively. The 
wave vector for each plot is in units of 20./ s, where s = speed of sound at the trap centre. Our QHR lattice is specified by 
fi/20 = 0.7. Triangle size reflects error bars, b) Extraction of the effective wavelength (A) for both branches of an azimuthally 
symmetric mode from the phase of the density fluctuation (p) of the bulk mode. In all phase plots we interpolate the bulk 
variations over the vortex cores and the phase changes from to 2ir as we move from dark to light. 




FIG. 3: a) Cross section (along z-axis) of the effective mean 
field potential (^mu^r 2 + 2grA r |$ | 2 — p) experienced by a 
particle in the trap b) The phase and c) amplitude of the 
Bogoliubov quasi-particle state (it) for this excitation. Rl = 
(r) (dashed circle) for a state in the Landau state fto.m (r, (f>) 
where m is chosen to match that of u in region I. 



deGennes mean field theory [l3l | in a pancake 
trap geometry. The ground state in a rotating 
frame is obtained from the Gross-Pitaevski functional 



H = J* 



1, 



QL, 



when 



h± = — + ^muj'^r 2 is the single particle Hamil- 

tonian in the x — y plane. We take L z as the angular 
momentum operator defined in polar (r, 4>) cordinates as 
—ih-J^. The quasi-two dimensional coupling constant 

g-2D = 27rg3£)/a z where g^n — (4ivaH )/m and a z is 
the oscillator length along the z-axis: a z = \Jh/mu> z - 
Unless stated otherwise we express energy, length and 
time in oscillator units: hui±, ^Jh/moj± and l/u>± re- 
spectively. Writing the mean field fluctuation as <5$ = 
u(r)e~ luJt — v(r)*e luJt , a first order expansion around the 
groundstate $ Q yields the coupled equations for the Bo- 
goliubov amplitudes u(r) and v(r): 



(h ± - nL z + 2.g 2D A|$ | 2 - n) Vj 
(h ± -nL z + 2g 2D N\<P \ 2 - 



fi) Uj - gN<P 2 Vj = ujjUj(l) 



We solve these equations numerically. In Fig. we 
plot resulting dispersion curves, for parameters relevant 
to the recently measured 0, Q azimuthally symmetric 
{Pt 0) 0iQ modes in the three physical regimes described 
above. Here p is the radial order of the excitation. As 
shown in (Fig. E>), we can obtain the effective wave 
length (A) for the excitation from a phase plot of the 
complex density fluctuation amplitude for the jth exci- 
tation pj = <&* u.j — $ Vj. Apart from quantum fluctu- 
ations, the real part of p has the same meaning as the 



hydrodynamic density fluctuation |13| . Consider the two 
scales labeled lot and toj for the branches in each plot. 
In accordance with hydrodynamic predictions [3 , the re- 
duction in C\ and Ci as the QHR is approached causes 
the inertial and Tkachenko waves to occur within increas- 
ingly smaller ranges in 10 around u> — and uj — 2. The 
triangular data point in Fig. represents recent JILA 
data la, la, in excellent agreement with this and related 
work 

To introduce the surface wave excitations - the BEC 
edge states - we employ the analogy between the BEC 
and an electron confined in a magnetic field [11]. We 
will see that the quasi-particle states (u) responsible for 
these surface waves map onto single particle states of 
the Landau level picture. First we rewrite the Hamilto- 
nian as h± — £IL Z = Hl + h(ujj_ — Q)L Z where Hl = 
(-ifiVi — mut±z x r) /(2m). Hl is equivalent to the 
Hamiltonian of an electron of charge e and mass m mov- 
ing in the x — y plane under the influence of a mag- 
netic field Bz with vector potential A — \Bi x r. We 
then make the identification eB/2m = lo±. If we de- 
note the eigenfunctions of Hl as h n ^ m with eigenvalues 
e„ jm = hu>±(2n+ 1), then the integers n and m represent 
the Landau level and the angular momentum quantum 
numbers respectively. The eigenvalue of — D,L Z may 
be written as uj n ^ n — h ((lj± + Q)n + (lo±_ — 0)rn + u>±). 
Note that the positive term — fi) m unilaterally lifts 
the degeneracy of the Landau levels and represents only a 
small perturbation in the rapidly rotating limit Q w lj±_ . 
To explore the analogy, we make two observations: First, 
the high angular momentum of these states ensures that 
|/c|£ 3> 1 such that the excitation is particle- like; that is 
that |i>| <C |w| implying that the properties of the exci- 
tation are dictated by the character of u. We may sepa- 
rately solve for u in the regions r > R and r < R (labeled 
I and II in Fig-EJ) from Eq.Q](top), ignoring the coupling. 
Second, the large momentum of these states confines \u\ 
mostly to region I, near the condensate surface. There, 
the effective potential is harmonic because the interaction 
term 2g2DN\<& \ 2 is negligible. Thus we find that the ex- 
citation coincides with the solutions of Hl represented by 
h n ,m- To visualize the edge states, we show in Figs. [31.:>,c 
a full (no approximation) numerical solution of |m| for an 
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FIG. 4: We plot the N p , l x and 8 denned in the text for different ground states. The labels Tm, Im, Sm approximately mark 
the locations of Tkachenko, inertial and surface modes respectively. The depletion is multiplied by a factor of 10000/iV. Note 
the Gaussian shape of the cloud in the QHR. 



example. We see that the amplitude is radially peaked 
near = \/2ml, with a spread of order I — l/Vfi. 
This is exactly as would be expected for the equivalent 
Landau state defined by h 0>m = ^i^r™ e lm<t 'e~ r ' '> /41 ' 
[To). We stress that these edge states cannot be ade- 
quately accounted for in any long wavelength treatment. 
Further, they cannot be neglected as they are bound to 
occur throughout the spectrum, and in all three regimes. 

Having identified the edge states, we turn to the spec- 
trum evolution across the three regimes. Our results are 
summarized in Fig. 01 We calculate the total excitation 
number N p which is related to the density of states \ 
by N p (uj) = f^xi^')^'. The second order contribution 
to the angular momentum l z = J(u*L z Ui — v* L z Vi)dr 
is due to the density and velocity field fluctuations, and 
Si = J\vi(r)\ 2 dr is the condensate depletion of the ith ex- 
citation at T — 0. By correlating these plots with MPEG 
movies of each excitation, we can assign a mode character 
to each individual spectral point. 

Consider the spectrum in the stiff regime. Particle- like 
excitations occur in the l z plots as families of points that 
organize into distinct lines with positive slope (labeled by 
Sm). By contrast, the more scattered points correspond 
to Tkachenko waves (Tm) in the interval < lu < 1 and 
to inertial modes (Im) in the interval u> > 2. As we move 
to the soft and then QHR, the l z spectrum is increasingly 
dominated by particle-like excitations mainly for two rea- 
sons. First, the associated drastic decrease in G\ and C2 
with increasing £1 affects the mode energies and pushes 
the Tkachenko and inertial modes into smaller intervals 
of lu around to — and ui = 2 respectively (see also 
Fig-EK)- For the same reason, these modes experience a 
reduction in their speeds and are shifted towards smaller 




FIG. 5: a) Amplitude and b) phase of p for an inertial wave 
of long wavelength, c) Amplitude and d) phase of p for an 
inertial wave of shorter wavelength. 



l z . Second, what is key is that as — * uj±, \k\£ — * 00 and 
inertial modes become increasingly organized and acquire 
particle- like character such that, deep in the QHR, they 
are well described by a superposition of closely spaced 
(single-particle) Landau level states h n ^ m . 

The progression from collective to single particle char- 
acter is well illustrated by the behavior of the ampli- 
tude and phase of p for £1 — > uo±_ as shown in Figs. |S] 
The mode in Figs. E^jb is consistent with a hydrody- 
namic approximation 0, 0, Q which stipulates a func- 
tional form p — f(r)e lp ^, where p is an integer. As J! 
increases and the lattice becomes softer, the approxima- 
tion breaks down for increasingly longer wavelengths and 
features signalling particle character begin to appear as 
interference structures on the scale of e (the 'islands' in 
Figs. Ek,d). For particle-like excitations, l z has a similar 
meaning to m for the h n , m states. We may thereby inter- 
pret the structure of the particle-like sub-spectrum from 
the expression for uj n ,m- If we identify the groundstate 
energy p, with the zero point energy (huj±) of w„_ m then 

the eigenvalue for h± — L z fi suggests bands centered at 
lu = 2n : n G {0, 1, 2, ..}. Physically, the large m of the 
surface states ensures that they occupy larger intervals 
of the lu axis around the Landau levels than inertial or 
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FIG. 6: a) Amplitude of p for a Tkachenko wave, b) and c) 
represent for degenerate bound Hartree-Fock states, where 
ip is the wavefunction. d) of a extended Hartree-Fock 
state. In all plots the overlayed white dots are the equilibrium 
vortex positions 

Tkachenko modes due to the term {lo±_ — f2) m. Thus, the 
whole spectrum acquires a band-like structure as seen in 
the N p plots of Fig. H 

Having identified the connection of the BEC surface 
waves to IQHE edge states, a casual glance at \p\ for a 
Tkachenko (Fig. and an inertial mode (Fig. [SJi) sug- 
gests an extension of the BEC analogy to other single par- 
ticle states of the FET: The Tkachenko waves, confined 
to the vortex cores, are analogues of the bound states 
and the inertial modes, which cover the entire array (see 
Fig-Eb); are analogues of the extended states. Conceptu- 
ally, the vortex cores play the role of charged impurities 
in the FET. However one important difference is that in 
a BEC, the Tkachenko modes are generally collective and 
not particle- like due to significant mean field dynamics. 
To explore the connection we solve the related Hartree- 
Fock problem defined by Eq.n mmu s the coupling term. 
This is equivalent to assuming a nearly incompressible 
lattice, i.e. to infinite values of C\ and Ci. In addition to 
surface excitations, we obtain both a group of degenerate 
single particle states bound to the vortex cores (Fig.|S|3,c) 
as well as extended states (Fig.|S|i) which cover the array. 
As the C\ and Ci acquire finite values, these states map 
back onto the Tkachenko and inertial modes of the full 
problem and have lower energies than their Hartree-Fock 
counterparts. 

Consider next the condensate depletion S. As we ap- 
proach the QHR, the density drops, |fc|£ — > oo and 
\k\e — > 0. The inertial modes become particle-like and 
contribute less to S, while the Tkachenko waves become 



more 'collective' and contribute more to S. Moreover, 
the total depletion (X^i) grows. This can be seen in 
Fig. 01 In concurrence with work predicting a quantum 
phase transition 0, ^| as fl — > u>±, we see in Fig. 
a significant increase in the depletion of the Tkachenko 
modes as fl — > u>±. Like phonon vibrations in a crystal 
lattice, the Tkachenko waves arise out of a broken con- 
tinuous symmetry, and are a path to the destruction of 
the array and, in this case, destruction of the condensate 
as thermal and/or quantum fluctuations grow. 

The IQHE framework for describing the excitations 
can be used to understand many properties of the ro- 
tating BEC. For example, vortex nucleation. A small 
increase in Q is equivalent in the electron picture to an in- 
crease in the magnetic field strength. As B increases, the 
radius of the classical orbital 'ring' is reduced to keep the 
flux enclosed by the degenerate electron orbit constant. 
In the BEC, the decrease in the edge state radius with 
increasing causes vortices of the precessing ring to be 
drawn towards and included into the interior of the con- 
densate. This evolution was detected in recent models in 
which vortices were observed to enter into the bulk of the 
condensate via a precessing ring ,14]. Quantitatively, we 
consider a perturbation of EqsQ]by fl — > fl + dtt starting 
from ~ 0. We find that the lowest edge state acquires 
a negative energy, falling below the ground state - a sig- 
nature of a thermodynamic instability. This mechanism 
is particularly relevant to the experimental method pio- 
neered at JILA 9] in which the condensate is nucleated 
from a rotating thermal cloud. 

In sum, the Landau level analogy provides an organiz- 
ing principle for vortex lattice investigations in rotating 
BECs. In addition to describing the particle-like exci- 
tations of the lattice which grow in number as £ — * oo, 
it inspires further analogies with single particle states of 
the FET. 
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